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Abstract. Gamma-ray bursts (GRBs) have long been considered as candidates of ultrahigh-energy cosmic rays (UHECRs). 
We investigate the signatures of CR proton acceleration in the GRBs by consistently taking into account their hadronic and 
electromagnetic interactions. We discuss the implications of our findings for high-energy gamma ray observations with the 
recently launched Fermi Gamma-ray Space Telescope. 
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INTRODUCTION 

Gamma-ray bursts are the most powerful explosions in the universe releasing about 10 51 ergs of energy in keV-MeV 
y-rays within tens of seconds. Emission of this non-thermal radiation via synchrotron and/or Compton scattering 
requires acceleration of electrons to very high energies, probably by a Fermi mechanism in relativistic shocks. Rapid 
variability time scales (t v ~0.001-1 s) observed in non-thermal y-ray light curves implies that the emission region is 
compact and moves with a large bulk Lorentz factor (T ^100-1000) towards the observer. 

In the context of the fireball shock model (see, e.g., Refs. JHI2l for reviews), successively ejected materials in the 
form of fireballs from a central engine (black hole or a magnetar) collide with each other, due to differences in their 
bulk Lorentz factors, forming relativistic forward and reverse shocks (internal shocks). A baryonic (assumed mostly 
protons) contamination/load in an otherwise pure e and y dominated fireball or the energy release from the central 
engine may change in time and lead to a varying bulk Lorentz factor. Thus a baryonic component is required in this 
model to explain rapid y-ray variability. 

Protons can be accelerated to UHE in the same shock which may be responsible for electron acceleration in the 
GRB jet yfl. The energy losses and size scale of the shocked region limit, however, the maximum proton energy. 
Synchrotron radiation and photohadronic (py) interactions are dominant energy loss channels for the UHE protons 
and the lost energy is converted to high-energy y-rays and neutrinos as sb i a nam. The acceleration of 
protons to UHE takes place on a longer time scale than the acceleration of electrons in the shocks. Therefore, a 
characteristic time delay is expected for proton-induced high-energy y-ray and neutrino emission as compared to the 
keV-MeV emission by electrons via synchro-Compton mechanism(s). We discuss these signatures of UHE cosmic ray 
acceleration in the GRB jet as currently being probed by the Large Area Telescope (LAT) and the Gamma-ray Burst 
Monitor (GBM) onboard the Fermi Gamma-ray Space Telescope. 



PROTON ACCELERATION AND ENERGY LOSSES 

The y-ray spectrum in the keV-MeV range is well-fitted by the Band function 1 12] which is defined as 

*(p\ (£/100keV) a exp(- e / eo ) ; e < (a-j3)e 

{ > \[(a-j3)(e /100keV)] a -' 3 exp(j3-a)(£/100keV)/ 3 ; e>(a-j3)£o, K ) 

where e is the photon energy measured in keV and srf is the normalization factor measured in units of pho- 
tons cm~ 2 s _1 keV -1 at 100 keV. An integration of the quantity eBUje) over a given energy band (which we assumed 
to be 1 keV - 1 GeV throughout our calculation) and time interval gives the fluence 5? (e.g., in units of ergs cm -2 ) 



TABLE 1. Derived quantities from Band function parameters: a = — 1, |3 = —2.2, 
Eq = 10 3 keV, si = 10~ 2 ph cm" 2 keV" 1 at 100 keV, and z = 2 for a 10 s long GRB 
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in that energy and time range. The total isotropic-equivalent energy release from a burst at redshift z and luminosity 
distance d L is S = A%d\5?l{\ +z). 

The proper number density of photons (e.g., in units of photons cm~ 3 keV -1 ) in the comoving fireball frame (all 
variables are primed in this frame) can be obtained from the Band function as 

where e' = e(l + z)/T and the size of the emitting region is R' m Yt v c/(\ +z). The energy density of photons in 

the comoving frame is UL = J , max e n Je )de . The energy density of baryons is assumed to be U' = UL/e e , where 

e e < 1, in the fast-cooling scenario. The magnetic energy density is U' B = B' 2 j%% — £bU' p , and the magnetic field is 

B' = ^j87t(£ B /£ e )Uy with e B < 1. We choose the typical Band function parameters a = — 1, j8 = -2.2, £o = 1 MeV, 

srf = 10~ 2 ph cm -2 keV -1 and z — 2 (corresponding to dt — 4.85 x 10 28 cm for the standard ACDM cosmology), and 
we have listed in Table [TJ different derived quantities for t v = 10~ 3 -1 s and for 1 keV < £ < 1 GeV. We assumed that 
r is the same as r m ; n which has been derived using the condition that the fireball becomes optically thin to 77 — > e ± 
pair production at 1 GeV (see, e.g., Refs. dHQ). 

The acceleration time scale for protons in the electric field induced by B' is proportional to the Larmor time (E' p /eB r ) 
and is given by 

'p.acc ~ m 2 c 4 ft' P ~ B\ ' K ' 

where B CI = m 2 p c 3 /eh = 1.488 x 10 20 G, B' = 10 4 B' A G, E' p = l0 9 E' p 9 GeV and <j) > 1 is the number of gyro-radii 
required to increase a particle energy by a factor of 2.7. The synchrotron cooling time for protons in the same magnetic 
field is 

Ar e m e cB' 2 E' p ^ B' 2 E' p /- {) 

The photohadronic cooling time scale, which does not depend on the magnetic field but on the observed keV-MeV 
y-ray source density in equation (O, is given by 

m 2 c 5 z 100 n' <e') r^E'e' lm„c 2 r oo 

^ 1= 2kL d£/J ^L de' r a PY (e' r )M)e^ca (5) 
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for single py scattering. We approximated a py cross-section of Ob = (T p y(e' r )f(e' r ) = 68 jib for e' r > £ t ' h following 
Ref. 01. Here e' r = e'(l - j3 p cos G)E p /m p c 2 is the photon energy evaluated in the proton's rest frame for the angle 9 
between the directions of the proton and target photon, and we take £ t ' h ~ 0.2 GeV as the threshold photon energy for 
pion production in the rest frame of the proton. 

We have plotted in Fig. [TJ different time scales for protons with the assumed Band function parameters, redshift 
and four different variability time scales (see Table [TJ. Two sets of three straight lines, with positive and negative 
slopes, correspond to acceleration and synchrotron cooling times respectively, for £b/£<, =10 (dashed lines), 1 (solid 
lines) and 0. 1 (dot-dashed lines) in Table [TJ The horizontal solid line in each panel correspond to the dynamic time 
f dyn ~ f vr/(l +z). The vertical thick-dashed line is the proton escape time t' esc = (3/2)(m 2 c 4 /ti)(B' /B CY )t'^ n /E' p in the 
Bohm diffusion limit, plotted for the proton energy at which fg SC = t' iyn . Note that for all four t v plotted in Fig. [TJ 
synchrotron cooling effectively determines the maximum proton energy if Eg/^e ~ !■ The photohadronic cooling 




FIGURE 1. The acceleration and cooling times (labeled in the lower-right panel) for protons for the GRB parameters assumed 
in Table [T] with different t v . The dashed, solid and dot-dashed curves in the sets labeled "acc" and "syn" are for Eg/Ee =0.1, 1 
and 10 respectively. We used the Band function in equation Q} in the 1 keV - 1 GeV observed photon energy range with derived 
parameters given in Tablefj] and <j> = 1 in equation l|3} for our calculation. The maximum proton energy can be read off from a plot 
when the "acc" curve, for a given B' field, is intersected by one of the cooling curves "py" or "syn" for the same B' field, or by the 
"dyn" curve. The top left panel for ty — 1 s and Bg/'Bg — 1 correspond to Ep : max — 

3 x 10 20 eV, the highest among all choices of the 

parameters and Eg/Ze- See main text for details. 



(smoothly-curved solid line labeled py) is important for highly variable GRBs as the internal photon density [see 
equation ©J increases to allow frequent interactions by high-energy protons. The protons may not escape the shocked 
region for the parameters assumed here. 

Proton-synchrotron radiation and photohadronic interactions (via 7t° decays) produce very high-energy photons. 
Photohadronic interactions also produce very high-energy electrons (via 7T+ decays). These particles initiate electro- 
magnetic cascades in the GRB fireball. We postpone the details of the cascade modeling and emission for future work. 
In the next section we concentrate on an analyic estimate of the proton-synchrotron radiation, which dominates for 
eg/fie ~ 1 in the GRB fireball. Note that, in this limit the Compton parameter Y = (— 1 + \/l +4£ e /£fl)/2 < 1 in the 
fast-cooling scenario and the Compton scattering of synchrotron photons to high-energy y-rays (SSC) may not be 
important flrfl . 



PROTON SYNCHROTRON RADIATION 

Assuming that synchrotron radiation limits proton acceleration, the maximum shock-accelerated proton energy as a 
function of time in the shocked fluid frame is 

0<t' <t' 

E'(t') = < / ^, 3 „ (6) 
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FIGURE 2. Maximum shock-accelerated energy (E pmaK ) and characteristic synchrotron cooling energy (E' pc ) for protons as 
functions of the comoving time. Note that the synchrotron cooling is important for t' >t' p0 where t' is denoted by the arrow. 



The maximum proton energy in the second case in equation (O is found from the condition t' p acc = t' p syn and the time 
to accelerate protons to this maximum energy are given by 

1 /2 

(9 tunic* B a\ ' 2xl0 9 
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respectively. Thus for t' > t' the maximum proton energy remains constant at E' p (see Fig. |2j. The synchrotron 
cooling, however, becomes significant with increasing time. The characteristic synchrotron cooling-break energy for 
protons from the condition t' p syn = t', and for t' dyn > t' > t' p is given by (see Fig. |2]l 

' 4 r e m e cB' L t B£t 2 

where t' 2 = f'/lOO s. Protons above this energy cool efficiently down to E'(t') within t' . As a result the proton 
spectrum can be approximated by a broken power law with the spectral index softened by unity for E' p > E' pc lfl7ll . 

The number density of particles behind the shock is 4n'r re i, where n' is the pre-shock number density in the GRB 
fireball. As the shock crosses the fireball, it sweeps-up particles and the number density of shocked particles increases 
°< t ' . We assume that a time-dependent fraction of the shock energy density An'm p c 2 T^ ei is channeled to the high energy 
protons. In the most optimistic scenario, the time-dependent proton number spectrum can be written as 

J ( E > t '\ ~ An'—c x / ( E 'p/ E P^ K ' E 'p-^ - E 'p < KA 1 ') ,, m 
n p[ t p ,t ) ^c, I ( Ep! jE pin J(E> p /E p!n J-*- 1 ; E' pc (t') < E' p < E> p>mla {t>) 

where C p is an overall normalization factor and E' min ~ T vs ,\m p c 2 is the minimum proton energy, and r re i ~ 1-10 is 
the relative Lorentz factors between two colliding shells. The snapshot of this spectrum is plotted in Fig.|3]at different 
t', and for arbitrarily fixed n' and C p . 

The synchrotron power and the typical synchrotron photon energy emitted by a proton of energy E' p are given by 

e' = (12) 
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FIGURE 3. Shock-accelerated proton spectrum at increasing comoving time (from bottom to top). Total number of protons 
swept-up by the shock increases linearly with time as in equation < | 1 01 >- At early times, t' < t' in equations $6$ and (8), the 

maximum proton energy increases linearly (dashed lines) and saturates to a constant value E' given in equation for t 1 > t' p 

(solid lines). As the synchrotron cooling becomes important for t' > t' , the spectrum shows a break at E' pc given in equation l|9} 

which decreases with time as oc . The spectrum is oc E'~ K below E' p c and oe E'~ K above E p c . 



respectively. For the characteristic proton synchrotron cooling energy in equation ©, the typical synchrotron photon 
energy from equation ( fT2l for t ! dyn >t'>t' p0 is 

243 h A R 3 1 22 
£ > oc ( t >) = , B . i« GeV. (13) 

°- cK ' 32 r 2 m 2 e m p c 4 B' 3 t' 2 B%t% 

The ratio of the synchrotron power to the typical photon energy roughly gives the spectral power. The total spectral 
emissivity (e.g., in units of ergs cm~ 3 s GeV -1 ) from all protons with energy E' pc from the distribution in 
equation ( fTOb is given by 

cm< tjs ^ P'iK.ci 1 ')) , , E , 8 r e m e m p <?B' 

The corresponding time-dependent energy flux (e.g., in units of ergs cm~ 2 s -1 ) of the proton synchrotron radiation, in 
analogy with electron synchrotron radiation in the slow-cooling regime 11711 . is 

J^,(0 - K^kS*') x { {£'/^Y {K - m ; C < e' < <c(0 (15) 

«c/^)" (K " 1)/2 (e'/<c)- K/2 ; <c(0 < e' < CM- 

Here e^ in and e' max are respectively the typical photon energies for the protons with energies E' p min and E p max . We have 
plotted this proton synchrotron radiation spectrum at different times t' > t' in Fig.O High energy photons e' > MeV, 
however, produce e ± pairs which initiate cascades. Note that the proton synchrotron radiation and the cascade radiation 
build-up gradually with time as the shock propagates through the GRB fireball. As a result, high-energy emission from 
proton synchrotron radiation would be delayed compared to the prompt keV-MeV emission from shock-accelerated 
primary electrons. 



DISCUSSION 



We have investigated different scenarios of fireball shock model for GRBs in which protons can be accelerated to 
ultrahigh energies and produce observable y-ray signatures. We found that — for the Band function parameters 
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FIGURE 4. Proton synchrotron radiation spectrum at three different times t' > e' oc . Note that the peak of the synchrotron 
radiation spectrum e' oc decreases with time as t'~ 2 as indicated by the arrows. The spectrum is °c gM^-i)/ 2 below e' oc and 
oc e'~ K / 2 above e' c> c , where k > 2 is the spectral index for shock-accelerated protons. 



adopted here, in the 1 keV - 1 GeV range, typical of a bright long duration GRB — synchrotron radiation is the 
dominating energy loss channel for UHE protons for the model parameters £s/£e ~ 1, the acceleration parameter 
<j) w 1 and for the variability time scale t v Z 10~ 3 s. An SSC component of high-energy y-rays, which are expected to 
form simultaneously with keV-Me V synchrotron photons, is suppressed in the case Eg/Ee ~ 1 as well. Photohadronic 
(py) interactions can dominate the energy losses by UHE protons if the GRB fireball is highly compact with ~1 ms 
scale variability or more target soft photons in the fireball than considered here, and for Eg/fie < 1 . In this case an 
SSC component is expected. We also found that proton synchrotron radiation and associated electromagnetic cascade 
emission, from synchrotron photons or/and pion-decay secondaries from p y interactions, can produce high-energy y- 
rays which are delayed compared to the prompt keV-MeV emission. The delay is due to the time required for protons 
to be accelerated to UHE when they suffer large energy losses. 
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